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Abstract 


We study the finite horizon Merton portfolio optimization problem in a general local-stochastic 
volatility setting. Using model coefficient expansion techniques, we derive approximations for 
the both the value function and the optimal investment strategy. We also analyze the ‘implied 
Sharpe ratio’ and derive a series approximation for this quantity. The zeroth-order approx¬ 
imation of the_ value function and optimal investment strategy correspond to those obtained 
by Merton ( 19691 ) when the risky asset follows a geometric Brownian motion. The first-order 
correction of the value function can, for general utility functions, be expressed as a differential 
operator acting on the zeroth-order term. For power utility functions, higher order terms can 
also be computed as a differential operator acting on the zeroth-order term. We give a rigorous 
accuracy bound for the higher order approximations in this case in pure stochastic volatility 
models. A number of examples are provided in order to demonstrate numerically the accuracy 
of our approximations. 


1 Introduction 


The continuous time portfolio optimization problem was first studied by Merton ( 19691 ). where he 
considers a market that contains a riskless bond, which grows at a fixed deterministic rate, and 
multiple risky assets, each of which is modeled as a geometric Brownian motion with constant 
drift and constant volatility. In this setting, Merton obtains an explicit expression for the value 
function and optimal investment strategy of an investor who wishes to maximize expected utility 
when the utility function has certain specific forms. However, much empirical evidence suggests 
that volatility is stochastic and is driven by both local and auxiliary factors, and so it is natural to 
ask how an investor would change his investment strategy in the presence of stochastic volatility. 
There have been a number of studies in this direction, a few of which, we now mention. Darius 


(2005) studi es t he finite horizon optimal investment problem in a CEV local volatility model. 
C hacko and Viceiral (2005j) examine the infinite-horizon optimal investment problem in a Heston- 
like stochastic volatility model. While both studies provide an explicit expression for an investor’s 
value function and optimal investment strategy, the results are specific to the models studied in 
these two papers and for power utility functions. 

Approximation methods, which have been extensively used for option pricing and related prob¬ 
lems, have been adapted for the portfolio selection problem, allowing for a wider class of volatility 
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models and utility functions. The Merton problem for power utilities under fast mean-reverting 
stochastic volatility was analyzed by asymptotic methods in (Fouque et al., 2 000 , Chapter 10), and 
the related partial hedging stochastic control problem in Jonsson a nd Sir car ( 20021)1 131) using asymp¬ 
totic analysis for the dual problem. More recently, Fouque et ah (2012) consider a general class of 
multiscale stochastic volatility models and general utility functions. Here, volatility is driven by 
one fast-varying and one slow-varying factor. The separation of time-scales allows the authors to 
obtain explicit approximations for the investor’s value function and optimal control, by combining 
singular and regular perturbation methods on the primal problem. These methods were previously 
developed to o btain explicit price approximations for various financial derivatives, as described in 
the book Fouque et ah ( 2011 ). 

Here we study the Merton problem in a general local-stochastic volatility (LSV) setting. The 
LSV setting encompasses both local volatility models (e.g., CEV and quadratic) and stochastic 
volatility models (e.g., Heston and Hull-White) as well as models that combine local and auxiliary 
factors of volatility (e.g., SABR and A-SABR). As explicit expressions for the value function and 
optimal investment strategy are not available in this very general setting, we focus on obtaining 
approximations for these quantities. Specifically, we will obtain an approximation for the solution 
of a nonlinear Hamilton-Jacobi-Bellman partial differential equation (HJB PDE) by expanding the 
PD E co efficients in a Taylo r seri es. The Taylor series expansion method was initially developed 
Pagliarani and Pascucci (. 2012 ) to solve linear pricing PDEs under local volatility models, and 


m 


is closely related to the classical parametrix method (see, for instance, Coriell i et ah 


(teoioi ) 


for 


applications in finance). The method was later extended in iLorig et all ( 2014bl ) to include more 
general polynomial expansions and to handle mu ltidimensional diffusion s. Additional l y, the tech¬ 
nique has been app lied to models with jumps; see Pagliarani et ah ( 20131 ). Lorig et ah ( 2014<j ) and 


Long et ah (j2014al ). We remark that the PDEs that arise in no-arbitrage pricing theory are linear, 


whereas the HJB PDE we consider here is fully nonlinear. 

The rest of the paper proceeds as follows. In Section 01 we introduce a general class of local- 
stochastic volatility models, define a representative investor’s value function and write the associ¬ 
ated HJB PDE. Section [3] presents the first order approximation and formulas for the principal LSV 
correction to the value function and the optimal investment strategy. Motivated by the notion of 
Black-Scholes implied volatility, we also develop the notion of implied Sharpe ratio, that provides 
a greater intuition about the resulting formulas, which are summarized in Section 13.71 We discuss 
higher order terms in Section 01 and show that power utilities are particularly amenable to ob¬ 
taining explicit formulas for further terms in the approximation. In Section 01 we provide explicit 
results for power utility. In particular, we derive rigorous error bounds for the value function in 
a stochastic volatility setting. In Section 01 we provide two numerical examples, illustrating the 
accuracy and versatility of our approximation method. Section [7| concludes. 


2 Merton Problem under Local-Stochastic Volatility 

We consider a local-stochastic volatility model for a risky asset S: 

^ Y t ) dt + a(St , Y t ) d B[ l) (1) 

d Y t = c(S t , Y t ) dt + p(S t , Y t ) d B [ 2) , 

where B ^ and B^ are standard Brownian motions under a probability measure P with correlation 
coefficient p £ (—1,1): E{d_B^ d B^} = pdf. The log price process X = log S is, by Ito’s formula, 
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described by the following : 


dX t = b(X t , Y t ) d t + a(X t , Y t ) d^ (1) (2) 

d Y t = c(X t , Y t ) d t + /3{X t ,Y t ) d B?\ 

where a(X t ,Y t ) = a(e Xt ,Yt), and similarly (//,c,/3) from (/2,c,/3), and we have defined 

The model coefficient functions (p, cr, c, /3) are smooth functions of (x, y ) and are such that the 
Markovian system (|2|) admits a unique strong solution. 


2.1 Utility Maximization and HJB Equation 

We denote by W the wealth process of an investor who invests 7p units of currency in S at time 
t and invests {Wt — ^ q) units of currency in a riskless money market account. For simplicity, we 
assume that the risk-free rate of interest is zero, and so the wealth process W satisfies 

d W t = ^ dS t = 7 T t n(X t , Y t ) d t + 7T t a(X t ,Y t ) d . 

The investor acts to maximize the expected utility of portfolio value, or wealth, at a fixed finite 
time horizon T: E{f7(Wr)}, where U : R+ —» R is a smooth, increasing and strictly concave utility 
function satisfying the “usual conditions” U'{ 0 + ) = 00 and U'{ 00 ) = 0. 

We define the investor’s value function V by 


V(t,x,y,w) := supE {U(W T ) \ X t = x,Y t = y,W t = w}, 

7rGll 


(3) 


where II is the set of admissible strategies vr, which are non-anticipating and satisfy 

e |' J n 2 a 2 (X t ,Y t )dt\ < 00 . 

and Wt > 0 a.s. 

We assume that V E C 1,2,2,2 ([0,T] x R x R x R + ). The Hamilton-Jacobi-Bellman partial 
differential equation (HJB-PDE) associated with the stochastic control problem ([3|) is 


¥r+A \ V + maxA n V = 0, 
at ) ttgR 


where the operators A and A n are given by 


V(T,x, y, w) = U(w), 


d 2 d 2 d 2 d d 

A = \a 2 {x, y )—j + pa(x, y)P(x, y)g^ + y)g^ + v)fa + c ( x > v)^ 


d 2 


A n = - 7 r a (x,y) — f + tt ( cr-(x,y) 


d 2 


dxdw 


+ pa(x,y)/3(x,y) 


8 2 


d 


ayBw +Mx ' y) aw 


We refer to the books Fleming and Soner ( 1993 ) and Pham ( 20091 ) for technical details. 
The optimal strategy 7r* = argmax^ AAV is given (in feedback form) by 


(4) 

(5) 


(a 2 (x,y)V xw + pa(x,y)/3(x,y)V yw + p(x,y)V w ) 
a 2 (x,y)V ww 


( 6 ) 
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where subscripts indicate partial derivatives. 

Inserting the optimal strategy it* into the HJB-PDE ((H) yields 


d_ 

at 


where 3ST(V) is a nonlinear term, which is given by 


+ A)V + N(V), V(T,x,y,w) = U(w), 


HV) = - 


(cr(x, y)V xw + p/3(x, y)V yw + X(x, y)V w ) 


2V W 


and we have introduced the Sharpe ratio 


\(x,y) = 


h{x,y) 

cr(x,y)' 


(7) 

( 8 ) 


2.2 Constant Parameter Merton Problem 


We review and introduce notation that will be used later for the constant parameter Merton prob¬ 
lem, that is, when y and a in ([TJ) are constant, and so therefore y and a are constant and the stock 
price S follows the geometric Brownian motion 

■= ydt + o <1B\ 1 ■* . 

St 

Then the Merton value function M ( t , w, A) for the investment problem for this stock, whose constant 
Sharpe ratio is A = y/a, is the unique smooth solution of the HJB PDE problem 

M 2 

M(T,w) = U(w ), (9) 

M-WW 

on t < T and w > 0. Smoothness of M given a smooth utility function U (as assumed above), as 
well as differentiability of M in A is easily established by the Legendre transform, which converts 
(|9|) into a linear constant coefficient parabolic PDE problem for the dual. Regularity results for the 
latter problem are standard. 

It is convenient to introduce the Merton risk tolerance function 

M 

R{t,w, A) := —(t,w;\), (10) 

M-WW 

and the operator notation 


V k :=(R(t,w,X)Y 


d k 

dw k ’ 


k = 1 , 2 ,-•• . 


( 11 ) 


We recall also the Vega-Gamma relationship taken from ( Fouque et ah . 20121 . Lemma 3.2): 
Lemma 2.1. The Merton value function M(t,w\ A) satisfies the “Vega-Gamma” relation 

d M 


d\ 


= -( T-t)AV 2 M , 


where V 2 is defined in dUD- 

Thus the derivative of the value function with respect to the Sharpe ratio (analogous to an 
option price’s derivative with respect to volatility, its Vega) is proportional to its negative “second 
derivative” DyM (which is analogous to the option price’s second derivative with respect to the 
stock price, its Gamma). This result will be used repeatedly in deriving the implied Sharpe ratio 
in Section 13.51 and the approximation to the optimal portfolio in Section 13.61 
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3 Coefficient Expansion &; First Order Approximation 

For general (c, /3, y, a, A} and U, there is no closed form solution of (0. The goal of this section is 
to find series approximations for the value function 

V = F (0) + F (1) + P (2) + • • • , 

and the optimal investment strategy 


7T — 7Tq + 7T| + 7T2 + 


using model coefficient (Taylor series) expansions. This approach is deve lo ped for the linear Eu¬ 


ropean option pricing problem in a general LSV setting in iLorig et al. ( 2014bh . where explicit 


approximations for option prices and implied volatilities are obtained by expanding the coefficients 
of the underlying diffusion as a Taylor series. Note that, here the HJB-PDE fl7|) is fully nonlinear. 
Our first order approximation formulas are summarized in Section [3.71 


3.1 Coefficient Polynomial Expansions 

We begin by fixing an point (x,y) E R 2 . For any function x that is analytic in a neighborhood of 
(x,y), we define the following family of functions indexed by a E [0,1]: 

OO 

y a (x, y) := Y an Xn(x , y), (12) 

71=0 

where 

n 

Xn(x,y) := Y Xn ~ k ’ k ' ( x ~ x) n ~ k (y - y) k , Xn-k,k ■= ^\ d x~ kd yX(x, y), 

and we note that xo = Xo,o = x(x,y) is a constant. Observe that x a | a =i is the Taylor series of % 
about the point (x,y). Here, a is an accounting parameter that will be used to identify successive 
terms of our approximation. 

In the PDE ([7]), we will replace each of the coefficient functions 

X € {At, c, ct 2 ,/3 2 , A 2 ,cj/3,/3A} 

by y“, for some a, (x,y), and then use the series expansion (1121) for x a ■ Another way of saying this 
is that we assume the coefficients are of the form 

x(s + a(x -x),y + a(y - y)j , 

whose exact Taylor series is given by GZD, and we are interested in the case when a = 1. 

Consider now the following family of HJB-PDE problems 

(^+A a )H a + N a (F a ) = 0, V a (T,x,y,w) = U(w), (13) 

where, for a E [0,1], A a and N a (-) are obtained from A and N(-) in ((SJ) and © by making the 
change 

{y, c, u 2 , /3 2 , A 2 ,cr/3, /3A} {y a , c a , (u 2 )“, (/? 2 )“, (A 2 )“, (<t/3)“, (/3A)“}. 
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The linear operator in the PDE (11311 can therefore be written as 


A a = Y^a n A n , 

n=0 


where we define 


•A-n ■— (W Jn 




d 2 


d 


d 


+ (\P 2 )n(x,y)-^ + b n (x,y)— + c n (x, y)— ,(14) 


dy 2 


dy' 


and the expansion of the nonlinear term is a more involved computation. 

We construct a series approximation for the function V a as a power series in a: 


V a (t,x,y,w) = Y,* n V {n \t,x,y,w). 

n =0 


(15) 


Note that the functions are not constrained to be polynomials in (x, y, w), and in general they 
will not be. Our approximate solution to (|7|), which is the problem of interest, will then follow by 
setting a = 1. 


3.2 Zeroth &c First Order Approximations 

We insert m into (TTTfl) and collect terms of like powers of a. At lowest order we obtain 


+ A 0 ) F (0) - 


<t q V^ + p(3 0 V y y Z + A 0 K) 


(o) 


do) 


= 0, V<®{T,x,y,w) = U{w), 


(16) 


where the linear operator Aq, found from (I14I1 . has constant coefficients: 


Q 2 Q2 q 2 & d 

Ao = f + P& oAjttw;—I- C °TT' 

ox- oxoy oy z ox oy 


(17) 


As a consequence, the solution of (fTBl) is independent of x and y: = V^°\t,w), and we 

have 


T/(°) 

y w 


V i 0) - = 0, V^(T,W) = U(w). 

Vww 


(18) 


We observe that (1181) is the same as the PDE problem ([9|) that arises when solving the Merton 
problem assuming the underlying stock has a constant drift jMj = y(x,y), diffusion coefficient 
Co = a(x,y) and so constant Sharpe ratio Aq = A (x,y) = y(x,y)/a(x,y). Therefore, we have 


w ) = M(t , w, Aq). 


The PDE (1181) can be solved either analytically (for certain utility functions U ), or numerically. 

Recall the definition of the risk tolerance function in (HOD and the operators T>k in (HID , where 
now we take in those formulas the Sharpe ratio Aq: 


R(t,w ; A 0 ) 


K 


( 0 ) 


y(°) 

y ww 


(t, w; A), 


D/,, 


A 0 )) fc 


d k 

dw k ’ 


k = 1,2,--- . 
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Proceeding to the order a terms in (fl3l) . we obtain 

+Al 0 ) y(1) + 5^2^ (1) +\ 2 0 V 1 vW + p/ 3 0 A 0 Pi|^« + m Pi^F ( 1) = -(iA 2 )^!^ 0 ). 

We can re-write this more compactly as 

(J- t + Aq + B 0 ) ^ (1) + Hi=0, F (1) (T, x, y, w) = 0, 

where the linear operator Bo and the source term H\ are given by 


(19) 


d d 

Bo = ^AqH 2 + A^Pi + pPoXqDi— + // 0 Z>i 


dx' 


Hi(t,x,y,w) = ( 5 A 2 )i(x,r/)Pil/ (0) (t )U ;). 


( 20 ) 

( 21 ) 


3.3 Transformation to Constant Coefficient PDEs 


Next, we apply a change of variable such that can be found by solving a linear PDE with 
constant coefficients. We begin with the following lemma. 

Lemma 3.1. Let V - ® be the solution of (|18p and let Alo and Bo be as given in m and (ED, 
respectively. Then satisfies the following PDE problem 

^+Al 0 + B 0 ) C (0) =0, V<®(T,w) = U(w). (22) 

Proof. This follows directly from observing that the nonlinear term in (1181) can be written 


v (o) 

v w 


y(°) 

V ww 


ci°J ) 


V^=V 2 vW, or 


Ci 0) ) ; 




vW = -V 1 V ( ®. (23) 


Therefore, from (1181) . we have 


d 


- + iA 2 P 2 + A 2 P 1 )c(°) = o, 


and (1221) follows from the fact that does not depend on (x, y). while Alo and the last two terms 
in the expression (1201) for Bo take derivatives in those variables. □ 

Next, it will be helpful to introduce the following change of variables. 

Definition 3.2. We define the co-ordinate z by the transformation 

z(t,w) = - log V^°\t,w) + \\l (T-t). (24) 

We have the following change of variables formula, as used also in ( Fouque et al. . 120121 . Section 
2.3.2). 

Lemma 3.3. For a smooth function V(t,x,y,w), define q(t,x,y,z ) by 

V{t,x,y,w) = q(t,x,y,z(t,w )). 
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Then we have 


— + Aq + j V — f + -A-q + Cq ) q, 


where the operator Cg is given by 


1,2 d 2 


d 2 


Co — 2^0^32 PPo^O “I" MO' 


d 2 


dydz dxdz 


(25) 


(26) 


Proof. We shall use the shorthand R^°\t,w) = R(t,w; Ao). From we have that 
and so, differentiating ([30]) . we find 

^ / y(o) \ ^ ^ 

V t = q t ~ ( ^5) +! A oU> V 1 V = q z , V 2 V = q zz -R$q z . 

Then, using the first expression in (1231) to write the PDE (fT 8 j) for E® as V"/ 0 ^ = |A qT> 2 V^°\ and 
differentiating this with respect to w gives 


v£ ] = \xl (rW) + A IrMr^vW. 

But from R^Vww = — eJ,°\ we have (R^) 2 Vwww = (Rw^ + 1)F^°\ and so 

v£ 0) = ±A 2 0 (rM + 1)VM - A 2 0 R^v£\ 

which gives that 


V, 


( 0 ) 


tw 


y(0) 

V W 


= -iA 2 0 (rW - 1 ). 


Therefore, we have 


+ 5 ^ 0^2 + AgX>i^ E — qt+ — 1) — ^Ag^ q z + 5 A 5 (q zz — r£ q£j + Ag q~, 


which establishes that 


More directly, we have 


d_ 

dt 


+ 7?AgT>2 + A5P1 ) v - ( — + ^Ag— ^ ) q. 


9 , 


dt 


dz 2 


(27) 


d 


d 


d 2 


d 2 \ 


p(3o\fDi — + )iqD\ — E = p/3 0 A 0 —— + ^0 I q, 

v dy dx J \ dydz dxdz J 

which, combined with (1271) . leads to (1251) . □ 

We define q ^ by V^°\t,w) = q(°\t, z(t,w)). Then the PDE (122|) for E^ is transformed to the 
(constant coefficient) backward heat equation for q(°\t,z): 

9 i 1^2 d 2 \ (g) _ g ' 

dt + 2X °d^r ~ 0 ’ 


qW(T,z) = U((U')- 1 (e- Z )), 


(28) 








but of course the transformation (|24j) depends on the solution V ® itself. Again, as does not 
depend on (x,y), while Ao and the last two terms in the expression (12611 for Co take derivatives in 
those variables, we can write 

+ Aq + Co^ = 0. (29) 

Now let q ^ be defined from by 

E (1) (t,z,y,u;) = q (1 \t,x,y,z(t,w)), (30) 

using the transformation (I24|) . Then, using Lemma 13.31 we see that the PDE (11911 for fA 1 ), which 
has (t, rc)-dependent coefficients through the dependence of "Bo in (l20ll on R(t, w; Ao), is transformed 
to the constant coefficient equation for q^: 

(j^+A 0 + C 0 )</ (1) + Qi = 0, qW(T,x,y,z) = 0. (31) 

The source term is found from Hi(t,x,y,w) = Qi(t, x, y, z(t, w)), where, from (l2lT) . we have 

Q!(t,x,y,z) = (±A 2 )i {x,y)qf\ (32) 

3.4 Explicit expression for V ^ 

In this section, we will show that V^\ solution of (I19H . can be written as a differential operator 
acting on First, we look at the PDE problem 

ftq + Q = 0, q(T,x,y,z) = 0, (33) 

where 9f is the constant coefficient linear operator 

A = — + Ao + Co- 
at 

We also suppose that the source term Q(t,x,y,z) is of the following special form: 

Q{t,x,y,z) = ^(T - t) n (x - x) k (y - y) l v(t,x,y,z), (34) 

k,l,n 

where the sum is finite and v is a solution of the homogeneous equation 


•Kv = 0. 

We first define (the commutator) Lx = [A, (x — x)I] by 

!K ((x — x)v) = (x — x)“Kv + Lxv , 


(35) 


(36) 


and so a direct calculation using the expressions (11711 and (l26ll for Ao and Co respectively shows 
that d d 

£y = (mo - 3*0 )I + + PVoPo-Qjj + Mo(37) 

where I is the identity operator. Similarly defining (the commutator) Ly = [A, (y — y)I] by 

ft ((y - y)v) = (y - y)ftv + L Y v, 
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leads to 


d d d 

&Y = cqI + 0q — + pcroPo-Q^ + Pfto ^0 -Q^ , (38) 

We next introduce the following operators indexed by s € [t, T]: 

M x {s) = (x-x)I+ (s-t)L x , My(s) = (y - y)I + {s -t)£y, (39) 

Then we have the following result by construction of these operators. 

Lemma 3.4. Recall that v solves the homogeneous equation ([35]) . Then 

9CM^(s)M(,-(s)n = 0, (40) 


for integers k,l. 

Proof. We first calculate 

TCM-xv = JAx'Kv + Lxv — Lxv + (s — t)TCLxv = (s — t)TCLxv, 

where we have used (1361) . But since J~C and L \ are constant coefficient operators which commute, 
we have tKLxv = Lx^tv = 0 using (|35j) . Therefore, given a solution v of the homogeneous 
equation, JVlxv also solves the homogeneous equation, namely (K (M^u) = 0. Iterating we have 
that = 0 for integers k. Similarly TOvlyV = 0 for integers l, and so the result (1401) follows. □ 

From this we can exploit the special structure of the source Q to obtain the following formula. 

Proposition 3.5. The solution to (1331) where the source Q is of the form (fMl) is given by 

q(t,x,y,z) = f (T - s) n M.x(s)M l Y (s)v(t,x,y,z)ds. (41) 

k,l,n ^ 

Proof. Due to the linearity of the problem, it suffices to consider a single term of the polynomial: 

Q(t, x, y, z) = (T- t) n (x - x) k {y - y) l v(t, x, y, z). 


Then, we check that the solution is given by 


q(t,x,y,z)= / (T — s) n M k x (s)'MY(s)v(t,x,y, z) ds 


by computing 


“Kq = —{T — f) n M^(t)My(f)i;(f, x, y, z) + J (T — s) n dfMY(s)My (s)v(t, x, y, z) ds 
= ~i T ~ t) n {x - x) k {y - y) l v(t, x, y, z) 

= -Q, 


using Lemma 13.41 for the second term. The formula (1411) in the general polynomial case follows, 
and clearly the zero terminal condition is satisfied by m- ° 

We can now solve for the first correction in the series expansion. 
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Proposition 3.6. The solution to m is given by 

q (1) {t,x,y,z ) = (T - t)X 0 A(t,x,y)qi°\t,z ) + \(T - tf X 0 B qf) (t, z), (42) 

where 


A(t, x, y) = Ai i0 [(x -x) + \(T - t)(p 0 - ±a%)] + A 0 ,i [(y - y) + ~{T - t)c 0 ] , 

B = A] ; oMo + XoapAoXo- (43) 

Proof. We observe that since q ^ satisfies the homogeneous PDE J~Cq^ = 0 from (1291) . so does qi° , 
which follows from differentiating the constant coefficient PDE for q(°>. Then applying Proposition 
13.51 with v = q^\ and n = 0, (k,l) € {(1,0), (0,1)} and substituting the definitions (l39l) for M.y 
and My leads to 

q (1 \t,x,y,z) = [(|A 2 )i i0 {(T — t)(x — x) + \{T -t) 2 L x ) 

+(AA 2 ) 0 ,i ((T - t)(y -y) + \{T - t) 2 L Y )] #(t, z). 

Finally, substituting for Lx and Ly from (|37l) and (|38|) and using that q ^ does not depend on 
(x,y) leads to (H2l) . □ 

In the original variables, this leads to 

EW (t , x, y, w) = (T — t)X 0 A(t , x, y) V 1 V<® (t,w) + \(T- t) 2 X 0 B V 2 V^ ( t,w ). (44) 


3.5 Implied Sharpe Ratio 

In an analogy to implied volatility, for a fixed maturity T and utility function U, one can define 
the Merton implied Sharpe rahic[{] corresponding to value function M(t,w ; A) of Section [2.21 as the 
unique positive solution A a (t,x,y,w) of 

V a (t, x, y, w) = M(t, w; A“). (45) 

The existence and uniqueness of the implied Sharpe ratio follows from the fact that (i) the function 
M satisfies M(t,w ) > U(w), since an investor with initial wealth w can always obtain a terminal 
utility U(w ) by investing all of his money in the riskless bank account, and (ii) the function M is 
strictly increasing in A. Since a higher implied Sharpe ratio is indicative of a better investment op¬ 
portunity, we are interested to know how local stochastic volatility model parameters {c, j3, //, <r, p} 
affect the implied Sharpe ratio. 

Using our first order approximation V a ~ + aVA\ we look for a corresponding series 

approximation of the implied Sharpe ratio as 

A a = A^ + aAW + • • • . 


Then, expanding 

M(t, w, A) = M(t, w ; A^°)) + aA^M\(t, w ; A^) + • • • , 
and comparing with the expansion 

V a {t,x,y,w) = M(t,w;X 0 ) + aV l - 1 \t, x, y, w) H- 

lr The authors thank Jean-Pierre Fouque for a number of fruitful discussions, from which the concept of the Merton 
implied Sharpe ratio arose. 
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yields A^ 0 ) = Aq and 


(46) 


A (i) = V (1 \t,x,y,w) 

M x (t,w; A(°)) ' 

Next, from Lemma EU we have 
M x (t,w, A®) = -(T-t)\ 0 V 2 M(t,w,\W) = -(T-t)X 0 V 2 V^(t,w) = (T - t)A 0 PiP (0) (L «;). 
Then, using the formula (14411 for V W in (1461) gives 


A (1 \t,x,y,w) = A(t,x,y) + \{T - t)B 


Vjv(°\t,w) 
V[ V(°){t,w)' 


By computing 


V 2 V(°\t,w) 
V 1 V(°)(t,w) 


Rw(t, w; A 0 ) - 1, 


we have 

A“ « A (0) + aA (1) = A 0 + a [A(t, x, y) + \(T - t)B ( R w (t , w, A 0 ) - 1)] , 
where R is the Merton risk tolerance function 


R(t, w, A) = 


M w (t,w; A) 
Alyjyj ( t , W, A) 


(47) 


3.6 Optimal Portfolio 

From Q, we have that the optimal strategy 7r“’* is given by 

a,* = v a (x,y)V£ _ pP a (x,y)Vy W _ V^_ 
(a a ) 2 (x,y)V“ w cr a (x,y)V* w V“ w 


(48) 


It is convenient in deriving a compact form for our portfolio approximation to write our first order 
approximation to the value function as the Merton value function evaluated at the first order series 
(1171) for the Sharpe ratio 

V a (t, x, y,w) ~V ( t , x, t, w ) := M(t, w, Aq + aA^^(f, x, y,w )). 


Then our approximate first order policy will be to substitute V for V a in (|48f) . 
We have 


V w (t,x,y,w ) 


V ww (t,x,y, w) 

V yw (t,x,y, w) 
V xw (t,x,y,w) 


M w (t,w; A 0 + aA (1 ^(t, x,y, w)^ + aM x (t,w; A 0 )^(T — t)BR ww (t) W] Ao ) + 0(a 2 ), 

M ww (t,w; A 0 + aA {1 \t,x,y,w)^j + a\(T - t)B {R ww {t,w, X 0 )M x (t,w, Ao)) w + 0(a 2 ), 

a\ 0t iM\ w (t, w ; A 0 ) + 0(a 2 ), 
a\i }0 M Xw (t, w; A 0 ) + 0(a 2 ), 


where 0(a 2 ) denotes series terms in powers of a 2 and higher. 

Let us compute 

V w _ M w (t,w\X 0 + aA^(t,x,y,w)) + aM x (t,w, X 0 )±(T - t)BR ww (t,w, X 0 ) 
Vww M ww (t,w; A 0 + aAA)(t,x,y,w)) + a\(T - t)B {R ww (t, w] X 0 )M x (t, w; X 0 )) w 
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= R (t, w; A 0 + aA (1) (t, x, y, w) \ - - t)BR ww + a\(T - t)B^-(R ww M\) w + 0(a 2 ). 

\ / 1V1 WW 1V1 W w 

= R (t,w; A 0 + aA^ (t, x, y, w)^j + a\(T - t) 2 BX 0 R 2 {R ww + RR WW w + (R w ~ 1 )Rww) + 0(a 2 ). 


Here we have used the following identities satisfied by the Merton value function M(t, w; A) and its 
risk tolerance function R(t,w; A): 


M x 


M„ 


M w 

M ww 


Mx 


M w 
M 2 

WVJ 


M Xw 


— (T — t)XR 2 , 

(' T-t)XR 3 , 

(T - t)XR 2 (R w - 1), 


where (|49|) comes from Lemma I2.lt ([50]) comes from multiplying (f49l) by 
expression (EH), we also use A! — (A ''w l)A/ly,. 

Additionally, we compute 


(49) 

(50) 

(51) 

—R] and in the last 


V, 


yw 


Vww 

Vxw 

Vww 


— “Aop 
— «Ai t o 


M\ w 

M ww 

M\ w 

M ww 


(t,w; A 0 ) 
(t,w, A 0 ) 


Therefore we have 


aA 0 ,i(T - t)X 0 R{t,w ; A 0 ) (. R w {t,w\ A 0 ) - 1) + 0(a 2 ), 
aXi fi (T - t)X 0 R{t, w, A 0 ) {R w {t, w\ A 0 ) - 1) + 0(a 2 ). 


7 r 


a,* 


ix a {x,y) 


(cr a ) 2 (x,y) 

+ a(T — t)X()R{R w — 1) 


| R (t,w, A 0 + aA (1) (t,z, y,w)) + ^a(T - t) 2 BX 0 R 2 {RR WW w + {R + R w ~ 1)#™™)} 
p(3 a (x,y) 


cr a (x,y) 


■Ao,i + Ago 


where R without an argument denotes R(t,w, Xq). One could substitute the first two terms of the 
polynomial expansion of the coefficients, but since they are assumed known, there is no loss in 
accuracy in using the full expressions. The first order approximate optimal strategy is written in 
terms of the risk tolerance function and its derivatives. 


3.7 Summary 

We collect here the expressions for our first order approximation formulas, which follow from the 
prior calculations and setting the accounting parameter a = 1. 

• Our first order approximation to the value function V(t,x,y,w) in Q, solution of the PDE 
problem © is given by V(t,x,y,w ) ~ V(t,x,y,w), where 

V ( t , x, y, w) = H (0) (t, w) + (t, x, y, w) 

= M(t,w; A 0 ) + ((T - t)X 0 A(t,x,y)T>i + \(T - t) 2 X 0 B V\) M(t,w, A 0 ), 

and A and B are given in (14311 . 

• The implied Sharpe ratio A = A (t,x,y,w) defined by V(t,x,y,w) = M(t,w, A) is approxi¬ 
mated to first order by A ~ A, where 

A(t, x , y, w ) = A(°) + A (1 ( 

= A 0 + A(t, x, y) + \(T - t)B ( R w (t , w\ A 0 ) - 1) • (52) 
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• Our first order approximation to the optimal strategy ir*(t,x,y,w) in ((6|) is given by n* « 7 r, 
where 

7 t(t,x,y,w) = ^ {-R ^0 + A (1) (f,z,y,ra;)) 

+ \{T — t) 2 B XqR 2 (RR WWW + (R + R w — 1)-Ruiuj)| 

+ (T - t)X 0 R(R w -1)( P/3 . (x, ^ A 0 ,i + Ai, 0 ) • (53) 

\v(x,y) ) 


This formula has principle term that is the classical Merton strategy --^R, but here is up¬ 
dated to account for LSV by using the current p(x, y) and a(x, y) values, and with the implied 
Sharpe ratio in the risk tolerance function. The other terms contain effects of correlation p, 
volatility of volatility (3, higher Taylor expansion terms of the stochastic Sharpe ratio, and 
higher derivatives of the the risk tolerance function with respect to wealth. Even for a utility 
function where there is no explicit solution for the constant parameter Merton value function 
M, the risk tolera nce is easily computed by numerically solving Black’s equation, as detailed 
in ( Fouaue et ah . 20121 . Section 6.2). 


4 Higher Order Terms 


Having obtained PDEs for and we examine the higher order terms. An exercise in 

accounting shows that for all n > 1 the function V^ n \t, x,y, w) satisfies a linear PDE of the form 


+ Aq + 23 0 ) V (n) + H n = 0, (T, x, y, w) = 0, 


(54) 


where the source term H n depends only on V^ (k < n — 1). To see this, observe that the rath-order 
PDE involves three types of terms 

1 


0(a n ): V t (n \ A k V^~ k \ ( k<n ), £ XjV^V^ 


V a 


, (55) 


where, in the last term, x is a place holder for one of the coefficient functions appearing in Af a , 
the symbols (a,<5) are place holders for (x,y) or null (meaning just a single derivative in w ), and 

— ) is the mth order term in the Taylor series expansion of ( t 4- ) about the point a = 0, i.e., 

V WW J JYI \ * VJVJ J 

k 


1 


V a , 


1 




( 0 ) 


+E' 


1 


V 1 


ww / k 


1 


E 


(- 1 )' 


Kw)k ~,(viS)>+"> L 


e n 1 

£lk,m j — 1 


fe) 

WW 


, (56) 


where Ik,m is given by 


Ik,m — {i — 


5 •> °m ) 


GN = (57) 

3 =1 

The terms in (|54ll that involve V^ are precisely those terms that appear in (^ + Aq + r Bo)V^ n \ 
The terms that do not involve are grouped into the source term H n . We provide here an 
explicit expression for the second order source term H 2 , which appears in the 0(a 2 ) PDE: 

= -a - <«i+®or«") • %+ A,v m - aw 


y(0) 

v WW 


y(°) 

V WW 


y(°) 

V WW 
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1,2N (vftotf)) , (vi 1) )(4 ) ) , (V^XT^) 


-2aA 2 x 


y(°) 

v WW 


- (p/3A)o- 


y(°) 

v WW 


- (p/3A)i 


y(°) 

v WW 


/I 2 p2\ \ v yw ) \Vw ) \ * xw ) \Vw )\Vxw ) 

- ( 2 P P JO —jn\ -PO-17775-Pi- 


y(°) 

v WW 


y(°) 

v WW 


y(°) 

v WW 


/ \ V XW )\Vyw ) ( \ 2\ \ v xw ) 

- (p°P)° —^- ( 2 ct )o_ ^r- 

*WW *ww 

Higher-order sources terms can be obtained systematically using a computer algebra program such 
as Wolfram Mathematica. 

Now let g(”) be defined from V^ by 

V^it.x^y.w) = q (n) (t,x,y,z(t.,w)), 


using the transformation (|24fl . Then, using Lemma 13.31 we see that the PDE (1541) for which 
has (t, rc)-dependent coefficients through the dependence of "Bo hr (T2Qh on R(t, w\ Ao), is transformed 
to the constant coefficient equation for q ( n ); 

+ ^0 + q^P + Qn = 0, q^ n \T,x,y, z ) = 0. (58) 

The source term is found from H n (t , x, y, w) = Q n (t , x, y, z(t , w)). 

We must establish that, for every n > 1 there exists a function Q n such that Q n (t, x,y, z(t,w)) = 
H n (t,x,y,w). From (1551) we see that the source term H n contains two types of terms, the first of 
which is AkV^-^ (1 < k < n ). Since Ak acts only on (x,y), we have that AkV^ n ~ k ^ = Akq^ n ~ k i 
The second sort of term appearing in (1551) are those of the form 


T/Wt/(0 m 1 _ r T / 

E En R 


j+k-\-l-\-m=n 

where we have used ([56]) . 
Next, using 


V, 


WW p= 1 


y( k )y( l ) o ( fe ) o (0 

v aw v 'yw HotzH'fZ 


y(°) 

v WW 




we see that (1591) can be written as 


l zG/fc p j 1 


L -*- y(°) 

— 1 V WW 


k, l, m < n — 1, 


-ctil , (<?i 0) + <tz))<tz 


y(°) 0 (°1 

y ww Hz 


+ 


(?i 


(0)^2 


(59) 


JQJO m ( p ( Jii) (JP) 4. 

qaz q^z S^ ( -1 \p I TT I ~ < lzz , Wz + qzz )qz 

x i (o) 2-X iJ Ev 11 (0) + , (1 

j+k+l+m.=n ( i z p =1 \i&Ik,p j =1 V l? 2 W* 


S °)\2 


(60) 


where ( l,k,m < n — 1). We have therefore established that, for every n > 1, the source term 
H n (t,x,y,w), which is composed of products and quotients of derivatives of V^ k \t,x,y,w) (k < 
n — 1), can be written be written as a function Q n , which is composed of products and quotients 
of derivatives of q( k \t,x,y,z ) (k < n — 1). 

In Proposition 13.61 we saw that qP\ the first-order transformed value function, can be expressed 
as a differential operator acting on q(°\ specifically X 1; = where 

A = [(T - t){\ A 2 )!(x, y)I + \{T - tf ((|A 2 ) li0 £x + (|A 2 ) 0 ,i£y)] (61) 
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We will show that, for certain utility functions U, each of the higher order terms q( n ' ) (n > 2) can 
also be written as a differential operator acting on q . From Proposition 13.51 we know that if the 
source Q n in the nth order PDE (15811 is of the form (1341) . then this will be the case. From (1321) . we 
see that 

Qi = Qi<?o where Qi = (^X 2 ) 1 (x,y)—, (62) 

Unfortunately, this is not always the case. 

To see this, we examine Q 2 , the source term in the PDE for q^ 2 \ which one can compute: 


//J 0 ) , „(ob (P « (1) \ 

<32 = (p 2 ) 2 9 f - - T 5 ) «3i + 6091 ) + Am 

V (?* ) Qx / 

( ( 1 ) \2 (!) (!) 

I i 1 ' 2 ! ^ qz ^ I o/l\2\ (1) , f m'. ^2 Qyz 

+ ( 2 A D—mi- nQz + \PPD 0—Tm - 

q z qz 

(i) (i) (i) (i) 

1 f a\\ (1) , /I 2 q 2\ q-yz qyz , <?2 qxz 

+ (pPtyiqyz + ( 2 P ft )o-7m-^ Po -7m— 

ql q z 

(1) (1) (1) (1) 

I (1) 1 / Q\ qxz qyz 2 \ qxz qxz 

+ PlTxz + (pvp)0 -Tm” + ( 2*7 )o-7m— ■ (63) 

qi qi 

From (1631) we see that Q 2 can be written as Q 2 = Q 2 £/o where 


0 /i A 2^ d ((qf ] + qfz)q [ z ] 

Q 2 = ( 2A %-|- 


q£z 


\ (S¥ Jyj(Qi + eoA)+^Ti 

+ + 2 (^ a2 )i^ £ i + (p/3A)o%^ 

qz qz 

+ {pf5X)l ^d~z Ll + ( ^ 2/32)o ^f r^d~ z Ll + 

rfi rp _ ( n D)\ fp 


(64) 


where Qi was given in (1621) . and £1 in (16T1) . 

In order to use Proposition 13.51 we must establish that coefficients of Q 2 are polynomials in 
(x,y,z). The complicating terms are those that contain derivatives of q (°) and q ^ divided by q% . 
Such terms are always polynomials in (x, y), but may not be polynomial in z. The following lemma 
provides conditions under which the differential operator Q n is guaranteed to have coefficients that 
are independent of z: 

Lemma 4.1. Suppose q(°\t,z) is of the form: 

q W (t,z) = a{t)e b ^ +zc ^. (65) 


Then, for every n > 1, the source term Q n appearing in PDE d58l) can he written as Q n = Q n q(°\ 
where the differential operator Q n has coefficients that are polynomial in ( x , y) and independent of 

z. 
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Proof. We will prove by induction on n that there exists a differential operator Q n whose coefficients 
are polynomial in (x,y), independent of z, and which satisfies Q n = Q n q^\ where Q n is the nth- 
order source term appearing in (1581) . We know from (1621) that such a Qi exists. We now assume 
such Qfc exist for all 1 < k < n — 1, and we show that Q n exists and has the required form. 

The existence Q& implies from Proposition 13.51 that there exists an operator L k such that 
q( k ) = £ k q(°). Moreover, since Q k is polynomial in ( x,y ) and independent of z it follows from (141 1) 
that L k has coefficients that are polynomial in (x, y) and independent of z. Now, we recall that 
Q n contains two types of terms: A k q l ' n ^ k> (1 < k < n) and terms of the form (j60l) . Let us first 
examine terms of the form A k q^ n ~ k ^ (1 < k < n). Note that 

A k q( n ~ k) = A k L n _ k q {0 \ 1 < k < n. 


The coefficients of A k are polynomial in (x, y ) and independent of z by construction. Hence, the 
coefficients of A k L n - k are also polynomial in ( x , y ) and independent of z. Now, let us examine the 
terms of the form (1601) . Using q^ = L k q we can express (EDI) as 


J2 * 

j-\-k-\-l+m=n 




( 0 ) 

ql 



+ q ( z°J)(dz£i j q (0) ) 




( 0 )' 


( 9 i 0) ) 


(0) 

Qz 


( 66 ) 


where l,k,m < n — 1. Since, by assumption, q^ is of the form (1651) . it follows that terms of the 
form (1661) are polynomials in (x, y) and independent of z. We have therefore established that Q n 
can be written as Q n q^ where Q n is a differential operator whose coefficients are polynomial in 
(x, y) and independent of z. □ 

We will see in the next section that, when U belongs to the power utility class, then is of 
the form (1651) . Thus, the nth-order term q^ n > can be written as a differential operator L n acting on 

q(0). 


5 Specific results for power utility 


In this section, we consider the case where the utility function U belongs to the power utility class 


Power utility: 


U{ w) 


w 1 7 
I-7’ 


7 > 0, 7 + 1- (67) 


For general LSV dynamics ((2|), we will obtain the second-order approximation for the value function 
u, optimal investment strategy n* and implied Sharpe ratio A. Then, we will establish error 
estimates for the approximate value function in a stochastic volatility setting. 


5.1 Value function 

To obtain the second order approximation to the value function V, we must first compute q(°\ q ^ 
and . With U given by (l67l) . we have U'{w) = u; -7 and [L/ r ] —1 (C) = ■ As found in lMertonl 

( 19691 ). we have 


w) = M(t , w, Ao) 


w 1 7 

--exp 

1-7 


1-7 

7 


A g(T-i)). 


( 68 ) 
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Then the transform variable z in (1241) is given by 

z(t,w) = rylogw + \ Ao(T - t), (69) 

and the solution of the heat equation PDE problem (1281) is 



Next, using (1421) and (l70l) . we compute 

q (1) (t,x,y,z) = (Mi.o (( t - 0 (®-*) + l( T -t) 2 (^o - ^o)) 

+ (^ A2 )o,i (( T_t )(y ~V) + U T ~ t ) 2 ( c o + ^PA)A 0 )) ^jq^°\t,z), 

from which we obtain 

V ( 1 \t,x,y,z ) = (i a2 ) 1>0 + \( T -tf (7P0 - 50-0)) ( 71 ) 

+ (iA 2 ) 01 ((T - t )(y-y) + l(T- tf (c 0 + 0 A 0 )) ) V^\t, z), 

where, as a reminder, (^A 2 )go and (^A 2 )o,i are given by 



Having obtained an explicit expressions for q and we can now compute g( 2 £ The second 
order source term Q 2 , given by (1631) . can be written as Q 2 = Q2Q0 where the operator Q2 is given 
by (1641) . From (1701) . we see that q is of the form (1651) . Thus, from Lemma 14.II we know that the 
coefficients of £>2 are polynomial in (x, y) and independent of 2 . Therefore, we can use Proposition 
13.51 to compute q (2 ^ = L 2 q 0 . The expression for q^ is quite long. As such, for the sake of brevity, 
we do not include it here. 

We can obtain from q^ using (l30l) and (l69l) . The same procedure can be used to compute 

higher-order terms: q^ (to > 3). Since the expressions for u ^ and higher-order terms are quite 
long, we do not present them here. However, in the numerical examples that follow, we do com¬ 
pute the second order approximation, and we will see that it provides a noticeably more accurate 
approximation of V than does the first order approximation. 


5.2 Optimal Strategy 


For power utility, the Merton risk tolerance function is especially simple: R(t,w, A) = w/ 7 , and 
it does not depend on t, T or the Sharpe ratio A. Therefore, the approximate first order optimal 
strategy in (l53l) is given by tt* rs tt, where 


n(t,x,y,w) 


p(x,y) 

a 2 (x,y) 


+ (T - t) A 0 



f pfi(x,y) 

V <r{x,y) 


Ao,i + A 



w 

7 


(72) 


which is also proportional to the current wealth level w as in the classical Merton strategy, but 
with proportion that varies with the model coefficients whose values move with the log stock price 
x and the volatility driving factor y. 
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We remark that it is also possible to compute the next order of the strategy approximation 7 r| in 
the case of power utility using the lengthy expression for . For the special case (x,y) = (x,y), 
we have 


7T 9 = w x 


(T — f) 2 (7 — 1) 
2y 3 


((7 - 1)(5A 2 ) 0 ,i(p/3A)i,o + 7(iA 2 )i )0 (50- 2 )i,o 


— (7 c i,o(5A 2 ) 0 ,i + ( 7 C 0 — (7 — l)(p/3A)o) A 2 )i 1 + 2p 0 (^X 2 )2fi - 7°o(5A 2 )2,o + (^A^i.oMi.o) 


+-— o 1/T~2\ —~ r(/ ocr / 5 )o(^A 2 )o,i (-2 (7C0 - (7 - l)(p/ 3 A) 0 ) (^A 2 ) 0 ,i + (~2p 0 + 7 °o) A 2 )i ;0 ) ) 

°i\2 a 7)o v ' 

5.3 Implied Sharpe ratio 

We now compute the first order approximation of the implied Sharpe ratio A, which was introduced 
in Section [3.51 From (|52|) . we have A ~ A, where 

A = A 0 + Ai i0 (x — x) + A 0 ,i(y — y) + ^(T-t)^A 0 ,i (c 0 + ^{pP A) 0 ) + Ai )0 (j^po - (^o)) )• ( 73 ) 

The second order correction A 2 is quite long, and we omit it for the sake of brevity. 

Observe that, for power utility, in which an explicit expression for the constant parameter 
Merton value function M is available, one can obtain an expression for the implied Sharpe ratio A 
by solving (1451) with M given by (1681) : 

A = / log (rR) (W )(r-i)' (74) 

This will be useful when we test the numerical accuracy of the Sharpe ratio approximation in two 
examples. 

5.4 Accuracy of the approximation for stochastic volatility models 

In this section, we establish the accuracy of 

n 

y(n) _ ^ V^ k \ 

k =0 

the nth-order approximation of the value function V, assuming stochastic volatility dynamics of 
the form 

dX t = (V(Yi) - ^ 2 (Y t )^j d t + cr(V t )dB t x , 

dY t = c{Y t )dt + p(Y t )dB Y , (75) 

d(B x ,B Y ) t =pdt, 

and a utility function U of the power utility class (1671) . 

Throughout this section, we will make the following assumption: 

Assumption 5.1. There exists a constant C > 0 such that the following holds: 

(i) Uniform elhpticity: 1/C < /3 2 < C. 

(ii) Regularity and boundedness : The coefficients c, pi3X, ft 2 and A 2 are C m+1 (IR,) and all derivatives 
up to order n are bounded by C. 

(iii) The risk aversion parameter in the utility function (l67l) satisfies 7 > 1. 
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Clearly, stochastic volatility dynamics (1751) are a special case of the more general local-stochastic 
volatility dynamics ([2]). As such, one can obtain a series approximation V^ of the value function 
V a (t,y,w), which is in this case independent of x, by solving the sequence of PDEs (I54j) . An 
alternative but equivalent approach is to linearize the full PDE (|13j) . and then perform a series 
approximation on the resulting linear PDE. This is the approach we follow he re. 

Assuming power utility (l67l) and dynamics given by (|75l) . Zariphopouloul ( 2001 ) shows that the 
function V a (t,y,w), solution of (fl3l) . is given by 


V a (t,w,y) = 


w 


1-7 


1-7 


(■ r(t,y)Y 


v = 


7 + (1 ~l)P 2 ' 


where the function ip a satisfies the Cauchy problem 

0 = (d t +A a )r, r(T,y) = 1, a €[0,1], 

and A a is a linear elliptic operator given by 

A a = (^ 2 ) a d 2 + (c a + i^(p/3A)“) d y + ^(iA 2 )“. 


(76) 


(77) 


(78) 


Let us denote A = A a | a= i and Y = 'iJA\a=i ■ 

Remark 5.2. Assumption 15.11 part (iii) guarantees that the last term in (178j) is strictly negative. 

Remark 5.3. The linearization transformation described above works only for one-factor pure 
stochastic volatility dynamics ([ 75 D . or complete market pure local volatility models, and only for 
power utility (1671) . For more general local-stochastic volatility dynamics ([2]) and utility functions 
U, one must work with nonlinear PDE ([71). 

We return now to m- Noting that A a can be written as 

A a = jr a n A n , A n = {\f 3 2 ) n d 2 y +(c n + ] —t(pp\) ^(iA 2 )„, ( 79 ) 

n= 0 V 7 J VI 

we seek a solution ip a to (|77l) of the form 

OO 

r = E ( 8 °) 

n =0 


Inserting (1 79 1) and (1801) into PDE (1771) and collecting terms of like powers of a we obtain the following 
sequence of nested PDEs: 


0(1): 0 = (d t + AcOV’o, Yo(T,y) 

n 

0(a n ) : 0 = (dt+A 0 )Yn + Yn(T,y) 


k= 1 


This sequence of nested PDEs has been solved explicitly in 
result here. 


Lorig et al. ( 2014bh 


= 1 , ( 81 ) 

= 0. (82) 

We present the 
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Theorem 5.4. Let -00 and ifn ( n > 1) satisfy (I8TT) and (1821) . respectively. Then, omitting y- 
dependence for simplicity, we have 


i’o(t) = exp f(T-t)~ —-(^A 2 ) 0 ) , 

V 77 z 

where the linear operator L n (t,T ) is given by 

n pT 


lfn{t) = L n (t,T)lf 0 (t), 


^ rl rl rl _ ^ ^ ^ 

£ n (f,T):=^/ dti J di 2 ■ ■ ■ / Sn(Mi)9i 2 (M2) ■ ■■Si k (t,t k ), 

fc=i ■' t ■' ti Jt k -i i n k 

with In k defined in dSZD and 

9i(t,t k ) = Ai(y(t,t k )), y(t,t k ) = y + (t k -t) ( c 0 T ~— ~(p/3 \)o ) + 2(A/3 2 ) 0 d y . 


Here, the notation Ai(ff(t,t k )) indicates that y is replaced by ^(t,t k ) in the coefficients of A{. 

Proof. See ( Lorig et al. . 2014bl . Theorem 7). □ 

Having obtained an explicit expression for if n (n > 0) we now define if n , the nth-order approxima¬ 
tion of if: 


ifn ■= y^k, with y = y. 


(83) 


k=0 


The accuracy of the series approximation if n is established in Lorig et al. ( 20131 1. 


Theorem 5.5. Let if be the solution of ([80]) with a = 1 and let if n be defined by (l83j) with ifi (i > 0) 
as given in Theorem\5-4\ Then, under Assumption^ 7J we have 


n -\-3 , 


sup I if(t,y)-if n (t,y)\ = 0 (t 2 ), 
y 


t := T — t. 


Proof. See (jLorig et all 120131 . Theorem 3.10). 


(84) 


□ 


Our task is now to translate the approximation if n and accuracy result for \if — if n \ into an approx¬ 
imation tA(”) a nd accuracy result for \V — V^\. Expanding V a , given by (1761) . in powers of a, we 
obtain 


V a = 


w 


1-7 


1-7 

re 1-7 


(if 


a\V 


1 — 7 U I-7 


k=1 \m=l 


<+( e —> (w) ( e p*. 

k,m J = 1 


= : E ak Vk, 


k=0 


where I k , m is defined in ([57]) . and 


n 1-7 1-7 k 1 / m \ 

, 7 (“) = E r (t >. r<«> = f—O v (t) = f— E ^ WO E n 7 • ( 85 ) 

fc=0 ^ ' m= 1 \*e/fc, m J=l / 

and y is set by y = y. The following theorem establishes the accuracy of V^- n \ the nth-order 
approximation of the value function V. 
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Theorem 5.6. Let (X,Y) have stochastic volatility dynamics (1751) and assume the utility function 
U is of the power utility class (ED. Then, under Assumption I5.il for a fixed w, the approximate 
value function , given by (f85D . satisfies 

sup | V(t,y,w) - V^ n \t,y,w)\ = O(r^), r:=T-t, 

y 

where r) is defined in (USD. 

Proof. Theorem 15.41 implies that ifo(t) = 0(1) as t —>• 0 and equation (fMD implies 

n _|_2 

sup i/j n (t,y) = 0(t~), n > 1. 

y 

It therefore follows from (1851) that V n satisfies 

n-\- 2 

sup 14 (t,y,w) = 0 (t~), n> 1 . 

y 


Therefore, we have 

sup | V(t,y,w) - V {n) (t,y,w)\ = 0 (t^), 
y 


as claimed. 


□ 


6 Examples 

In this section we provide two numerical examples, which illustrate the accuracy and versatility of 
the series approximations developed in this paper. Both are based on power utility, but the first 
order approximations described in Section 13.71 could be computed for utility f uncti ons outside of 
this class, for instance mixture of power utilities, introduced in Fouque et al. ( 2012 :. which allow 
for wealth-varying relative risk aversion. There the solution of the constant parameter Merton 
problem M is computed numerically, and LSV corrections in the formulas of Section 13.71 can be 
obtained by numerical differentiation. 


6.1 Stochastic volatility example 

In our first example, we consider a stochastic volatility model in which the coefficients (p, cr, c, j$) 
appearing in (J 2 D are given by 


t(v) = 


°(y) = 


i 


Vv' 


c(y) = K (0 - y), 


P(y) = 5y/y. 


( 86 ) 


Here, the constants (n,Q,5) must satisfy the usual Feller condition: 2nd > 6 2 . 

Assuming power utility (1671) . an explicit formula for the value function of the infinite horizon 
consumption problem is obtained in Chacko and Viceira ( 20051 ). For the terminal utility optimiza¬ 
tion problem that^we consider in this paper, an explicit formula for the value function V in Q is 
obtained in (Fouque et ah, 2012, Section 6.4): 


V(t,y,w) = 


w 


1-7 


1-7 


a r]A(T-t)y+r)B(T-t ) 


V = 


7 + (1 -7 )p 2 ' 


(87) 
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where 


B(t) = k9 


1 — (a_/a+)e at 
1 - (o_/a+) 


A(t) = a + 


1 - e -“* 

1 - (a_/a_)_)e _at ’ 



a± = 


p = -S 2 


—q ± \Jq 2 — 4pr 


2 p 


a = \J q 2 — Apr , 
'1-7' 




7 


PP~ K, 


r = 


1 ^ 1—7 

VI 


P~ 


An explicit formula for the optimal investment strategy 7 r* can be obtained by inserting (|87D into 
([ 6 ]). Likewise, an explicit formula for the implied Sharpe ratio A can be obtained by inserting (|87j) 
into ([Tip . 

The zeroth, first and second order approximations for V, it* and A can be obtained using the 
results of Section [5j Fixing y = y and using (f 68 |) . (1711) . (1721) and (l73j) . we obtain 


y(°) = 


* ~ 

7 r « 

Aq = 


U{w) exp > 

™py+™p5p 2 ^(T-t)y, 

7 7 7 

pVv, 


y(1) = (jk{ 0 - y) + (1 - , 


Ai = 


(T-t) 



where in the strategy we have also expanded the coefficients in Taylor series. The second-order 
terms are omitted for the sake of brevity. 

In Figure [U we plot as a function of er = 1/yfy the exact value function V, the exact optimal 
investment strategy 7 r* and the exact implied Sharpe ratio A. We also plot the zeroth, first and 
second-order approximations of these quantities. For all three quantities, we observe a close match 
between the exact function (u, n* and A) and the second order approximation. Figure [T| also 
includes a plot of the implied Sharpe ratio (both exact A and the second order approximation) as 
a function of the risk-aversion parameter 7 for three different time horizons. It is clear from the 
figure that the approximation is most accurate at the shortest time horizons, consistent with the 
accuracy result of Theorem 15.61 


6.2 Local volatility example 

We now consider a local volatility model in which the coefficients (y, a) appearing in (J2J) are given 
by 


y(x) = y, 


a(x) = 5e vx . 


( 88 ) 


Since Y plays no role in the dynamics of X, the coefficients c and f5 do not appear. 

Assuming power utility (1671) . an explicit formula for the value function V in this setting is 
obtained in Darius ( 20051 ) 


V ( t , x, w ) 


w 1-7 

1-7 


(/(t,e - 2 ^)) 7 


A(t) = e A+?? ( 2r?+1 )( T-t ) 


A_ - A + 

A +e 2< 7 2 (A+-A _)(T-i) 


277+I 
2 77 


f(t,d) = A(t)e B ^ d , (89) 

B (.t ) = 
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Figure 1: The value function (top left), optimal strategy (top right) and implied Sharpe ratio (bottom left) are 
plotted as a function of instantaneous level of volatility a = I/ffy assuming power utility (1671) and dynamics 
given by stochastic volatility model (1861) . In all three plots, the solid line corresponds to the exact function 
and the dotted, dashed and dot-dashed lines correspond to the zeroth, first and second-order approximations, 
respectively. The parameters used in these three plots are: T — t = 4.0, w = 1.0, k = 0.3, 9 = 0.2, <5 = 0.3, 
p = —0.75, p = 0.3 and 7 = 3.0. On the bottom right we fix the volatility cr = 0.3 and we plot the implied 
Sharpe ration A as a function of 7 for three different time horizons T — t = {1, 2,4} corresponding to black, 
blue and red, respectively. The solid lines are exact. The dot-dashed lines correspond to the second-order 
approximation. The parameters used in the bottom right plot are w = 1.0, k = 0.3, 9 = 0.2, <5 = 0.3, 
p = —0.75 and p = 0.3. 


m 


A + (l - e V(A + -A_)(T-^ 
l-(A+/A_)eV(A+-A_)(T-t)’ 


x Vat? 

* “ 2,7 


The optimal investment strategy ir* can be obtained by inserting (1891) into (O An explicit expression 
for the implied Sharpe ratio A can be obtained by inserting (l89|) into d74|) . 

The zeroth, first and second order approximations for V , n* and A can be obtained using the 
results of Section [5j Fixing x = x and using (l 68 l) . (1711) . (l72|) and d73j) . we obtain 


F (0) = U(w) exp - t)^j , W 1 ) = (1 7 ^ 2 e~ 2 ^ (// - 7 ±<5 2 e 2 ^) 


7 r 


w JL e -*nx _ w\- 2 nx 


'yS 2 


7 


5 2 


uo, 
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A 0 = ^e~ vx 


Al = —2^—y e 


rjx 


(m - 1 ¥' 


2 e 2?jx 


Once again, the second-order terms are omitted for the sake of brevity. 

In Figure [21 we plot as a function of a = 5e vx the value function V, the optimal investment 
strategy ir* and the implied Sharpe ratio A. We also plot the zeroth, first and second order 
approximations of these quantities. For all three quantities, we observe a close match between the 
exact functions (V, n* and A) and their second order approximation. Figure [2] also contain plots 
of the implied Sharpe ratio A (both exact A and the second order approximation A 2 ) as a function 
of the risk-aversion parameter 7 for three different time horizons. 




Figure 2: Value function (top left), optimal strategy (top right) and implied Sharpe ratio (bottom left) are 
plotted as a function of a = Se^ assuming power utility (1671) and dynamics given by local volatility model 
(1881) . In all three plots, the solid line corresponds to the exact function, and the dotted, dashed and dot-dashed 
lines correspond to the zeroth, first and second-order approximations, respectively. The parameters used in 
these three plots are: T — t = 5, w = 1.0, 7 = —0.8, S = 0.3, p = 0.3 and 7 = 3.00. On the bottom right 
we fix the volatility a = 0.25 and we plot the implied Sharpe ration A as a function of 7 for three different 
time horizons T — t = {1,3,5} corresponding to black, blue and red, respectively. The solid lines are exact. 
The dot-dashed lines are our second order approximation. The parameters used in the bottom right plot are: 
w = 1.0, 7 = —0.8, 5 = 0.3 and p = 0.3. 

7 Conclusion 

In this paper we consider the finite horizon utility maximization problem in a general LSV setting. 
Using polynomial expansion methods, we obtain an approximate solution for the value function and 
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optimal investment strategy. The zeroth-order approximation of the value function and optimal 
investment strategy correspond to those obtained by Merton (19691 ) when the risky asset follows a 
geometric Brownian motion. 

The first-order correction of the value function can always be expressed as a differential op¬ 
erator acting on the zeroth-order term. Higher-order corrections can always be expressed as a 
nonlinear transformation of a convolution with a Gaussian kernel. For certain utility functions, 
these convolutions can be expressed in closed-form as a differential operator acting on the zeroth- 
order term. Corrections to the zeroth-order optimal investment strategy can be obtained from the 
approximation of the value function. 

We also introduce in this paper the concept of an implied Sharpe ratio and derive an approxi¬ 
mation for this quantity. We obtain specific results for power utility and give a rigorous error bound 
for the value function in a stochastic volatility setting. Finally, we provide two numerical examples 
to illustrate the accuracy and versatility of our approach. The expansion techniques presented in 
this paper naturally lend themselves to other nonlinear stochastic control problems. Recent results 
for indifference pricing of options contracts have been developed in iLorigl ( 2014 ). 
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